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INTRODUCTION 


Understanding  of  the  dynamic  response  of  underwater  towed  systems  has  been  of  interest  for 
many  years.  In  response  to  this  need,  many  computer  programs  have  been  developed  over  the  past 
30  ye^  to  model  this  response;  these  coiiq>uter  programs  are  based  on  a  number  of  methods 
including  the  ctmtinuum,  finite  element,  finite  segment,  and  lumped  parameter  approaches.  Ref¬ 
erences  1  through  11  describe  the  underlying  methodologies  for  some  of  the  c(^s  in  existence 
today.  Although  the  calibration  and  validatimi  of  these  c^s  is  not  always  straight  forward,  it  is 
generally  believed  that  they  ci^ture  the  basic  dynamics  of  these  systems. 

Many  of  these  codes  generate  a  set  of  nonlinear  ^uations  of  motion  that  are  numerically 
integrated  to  determine  the  two  or  three  dimensional  motion  of  the  system  over  a  specified  interv^ 
of  time.  Because  these  codes  function  in  the  time-domain,  they  may  be  used  to  perform  "numerical 
experiments"  to  gain  insight  into  the  fundamental  dynamics  of  these  systems.  However,  this  is  at 
best  a  trial-and-error  process  because  the  equations  of  motion  general  and  solved  by  the  codes 
are  quite  omiplex,  and  in  themselves,  provide  few  useful  hints  about  the  ovoall  system  behavior. 

Alternatively,  some  of  the  codes'  break  down  dynamic  response  into  a  fewo*  number  of 
response  modes  of  the  system.  In  particular,  Refinence  1  discusses  die  classical  normal  mode 
reqxmse  of  the  system  based  on  a  numerical  linearization  of  the  nonlinear  equations  of  motion 
about  a  stationary  steady-state  equilibriunL  As  Reference  12  suggests,  dus  approach  offers  sig¬ 
nificant  insight  into  the  dynamical  response  of  the  system  and  into  how  ^s  response  is  affected  by 
various  system  parameters.  Also,  by  con^iaring  the  linear  and  nonlinear  responses  by  time  domain 
simulation,  ncmlinear  effects  can  be  identified  and  evaluated.  By  this  latter  procedure.  References 
1  and  5  recognized  the  nonsymmetric  response  of  a  mwed  sphere  under  tow  point  excitation  to  be 
a  result  of  system  nonlinearities. 


This  paper  presents  an  overview  of  a  procedure  that  combines  the  stricdy  numerical  finite 
segment  approach  msented  in  References  5  9  with  the  nmdinear  analytical  peatiobation  method 
of  multiple  scales” This  iqiproach  involves  a  two-step  process.  First,  die  fully-coupled  three- 
dimensional  nonlinear  equations  of  motion  for  a  finite  segment  model  of  a  towed  system  are 


decoriqxised  about  a  statitmary  steady-state  equilibrium  position  into  a  set  of  equations  in  nxidal 
(or  mnmal)  variables;  only  quadratic  nonlinearities  are  maintained  in  this  process.  These 
cmnplex-valued  equations  are  completely  uncoupled  in  their  linear  terms,  butin  general,  they  remain 
coupled  through  the  nonlinear  terms.  Inspection  of  these  equations  tells  which  modes  are  coupled 
at  dus  level  of  qiproximation.  Next,  the  method  of  multiple  scales  may  be  tqiplied  to  these  equations 
to  determine  the  effects  of  quadratic  ntmlinearities  on  the  response  of  single  modes  or  on  the 
interacti<His  between  naodes  of  the  system  under  tow  point  excitation.  This  report  outlines  the 
tqrplication  of  multiple  scales  to  the  response  of  a  angle  mode  only;  application  to  naode  coupling 
will  be  presented  in  a  future  report 


This  iqrproach  combines  the  benefits  of  both  numerical  and  analytical  methods.  Equations 
of  nootion  about  general  steady-state  ctmfigurations  can  be  generated  for  complex  systems,  and 
after  modal  dectmiposition,  analytical  procures  can  be  iqrplied  to  gain  insist  into  the  nxxlal 
response  of  the  system  about  those  ctmfigutations. 
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The  balance  of  this  report  is  divided  into  five  sections.  The  first  discusses  the  full  nonlinear 
equations  of  motion  of  a  towed  system  and  how  they  can  be  reduced  to  describe  motions  about 
stationary  steady-state  configurations;  terms  up  throu^  quadratic  nonlinearities  are  maintained  in 
the  prt^ss.  The  second  describes  the  modal  deramposition  and  solution  of  the  linearized  equations. 
The  third  discusses  a  set  of  reduced  nonlinear  equations  for  the  modal  variables  in  which  the  modes 
of  the  system  are  coupled  through  the  quadratic  nonlinearities.  The  fourth  section  presents  the 
methodology  for  the  application  of  the  metiiod  of  multiple  scales  to  determine  die  effects  of  the 
nonlinearities  on  the  response  of  single  underdamped  modes  of  the  system;  it  also  includes  results 
computed  for  two  modes  of  a  towed  sphere  with  excitations  both  in  and  normal  to  the  towed  end 
of  the  cable.  The  final  section  contains  a  short  summary. 


EQUATIONS  OF  MOTION 


NONLINEAR  EQUATIONS 

In  this  paper,  a  towed  cable  system  is  assumed  to  be  a  single  or  multiple  branched  cable 
system  with  towed  bodies.  The  system  is  assumed  to  be  towed  tom  a  single  point  The  cable 
baches  form  an  open-tree  system  having  no  closed  kinematic  chains.  Each  length  of  cable  may 
have  different  physical  properties,  and  the  towed  bodies  may  be  simple  spheres  or  more  sophisticated 
underwater  vehicles.  The  motion  of  the  system  tow  point  is  arbitrary.  Figure  1  depicts  a  single 
branched  system  with  a  towed  body. 

These  systems  maj^  be  modeled  by  segmenting  the  cable  into  a  series  of  discrete  rigid  links 
connected  end-to-end  using  fnctionless  spherical  joints.  In  general,  the  equations  of  motion  of 
these  so-called  finite-segment  models  noay  be  written  in  the  form 

(i  =  l,...,N;y=l,...,2N)  (1) 

0’  =  1,...,2N; k  =  1 . M)  (2) 

whm,  for  example,  the  jc,({  =  1, . .  .,N)  may  represent  the  orientation  angles  of  the  cable  links  and 
the  towed  vehicles  relative  to  some  reference  frame,  the  x,({  =N  + 1,  ...,2N)  represent  the  first 
derivatives  of  these  angles,  and  the  =  l,...,Af)  represent  extonal  inputs  such  as  motions  of  the 

system’s  tow  point  and  motions  of  the  towed  vehicles’  control  flaps.  There  are  many  different 
analytical  approaches  that  can  be  used  to  provide  the  equations  in  this  fonn;  the  methodology  used 
in  this  work  is  outlined  in  Reference  S. 

During  steady  forward  or  steady  (circular)  turning  motion,  and  in  the  absence  of  external 
disturbances,  a  towed  system  exhibits  stationary  steady-state  equilibrium  positions.  In  these  situ¬ 
ations,  Xid  =  1,  ...,N)  the  orientation  angles  of  the  links  and  towed  vehicles  relative  to  the  mean 
ship  frame  remain  constant.  In  this  case,  the  tight  sides  of  equations  (1)  and  (2)  are  zero,  and  the 
an^es  defining  the  equilibrium  shape  of  the  system  may  be  found  by  solving  the  N  nonlinear 
algebraic  equations 


2 
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0=0  (I J  =  k  =  1, (3) 

where  the  xj  and  ui  represent  the  equilibrium  values  of  the  orientation  angles  and  the  external 
inputs,  respectively. 

Numerically  solving  Equations  (3)  for  the  steady-state  equilibrium  orientaritm  angles  can  be 
difficult  unless  guesses  very  near  the  equilibrium  values  are  us^  to  initialize  the  solution  process. 
In  practice,  Aese  etjuations  can  be  repku^  by  a  series  of  smaller  sets  of  equations  that  can  be  solved 
sequentially,  allowing  for  a  much  more  rapid  and  effective  determination  of  the  i^uilihrium  angles. 
For  exanqile,  the  equSibnrium  equations  for  each  body  (cable  links  and  towed  vehicles)  in  the  system 
may  be  written  and  solved  as  separate  sets  of  equatimis  so  long  as  the  solution  process  begins  with 
equatitms  associated  with  bodies  at  the  ends  of  the  branches  and  then  proceeds  from  body  to  body 
towards  tire  system  tow  point  For  partially  submerged  systems  or  ^sterns  undergoing  circular 
modem,  dus  procedure  is  complicated  somewhat  In  tlrese  cases,  the  procedure  must  be  repeated 
until  convergence  for  the  over^  system  occurs;  however,  diis  iteradon  process  has  generally  been 
found  to  converge  very  rapidly. 


REDUCED  NONLINEAR  EQUATIONS  -  QUADRATIC  NONLINEARTTIES  ONLY 

To  describe  modems  of  the  system  that  result  feom  small  disturbances  to  the  state  variables 
Xiii  =  1,  ...,2N)  and  to  the  external  input  variables  Utife  =  1,  ...,Af ),  the  right  sides  of  the  nonlinear 
^uadons  of  modem  (1)  and  (2)  can  be  expanded  in  a  Taylem  Series  about  the  equilibrium  config- 
uradem  eleftned  by  the  soludon  of  Equadons  (3).  To  this  enel,  introduce  the  variables  y,  to  represent 
perturbadems  of  the  equilibrium  values  of  the  stare  variables,  and  introeluce  the  variables  to 
represent  perturbadems  of  the  equilibrium  values  of  the  external  input  variables  so  that 

jc,  =  Jc; + y,  (1  =  1,.. .,  2N)  (4) 


«*  =  u;+Vt  (*  =  l,...,Af).  (5) 

Subsdtudng  these  expressions  into  the  nemlinear  equadems  of  modem,  expaneling  in  a  Taylor  Series 
about  the  eqtulilnium  cemfiguradon,  and  exnitdng  terms  of  third  and  higher  order  in  the  poturbadons 
y,  and  results  in  equadons  of  the  form 


=  lA^yj + +^2  C^jy^ 

j  is  ^  ^  hpi  y»i* 


where 


A  -  —  «  —  _ 


(6) 

(7) 


c  D  =-^  E  - 


(8) 
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The  subscript  e  on  the  partial  derivatives  indicates  that  they  are  evaluated  at  the  equilibrium  con¬ 
figuration.  Because  the  first  N  equations  of  motion  (Equations  (1))  are  simply  definitions  of  the 
state  variables  jc,  (i  =N  + 1,..., 2W^),  many  of  the  Aj^,Bi(k,C^,/>iu,and£i(*  take  on  the  values  of  zero 
or  one  and  may  be  determined  by  inspection.  All  the  other  values  may  be  approximated  using  finite 
difi^nces;  in  this  work,  second-tmler  central  differences  are  us^  Also,  note  that  due  to  the 
equality  between  mixed  second  partial  derivatives,  the  C  and  D  arrays  possess  some  symmetry  so 
thatCj(,«  =  Ci^  (j,m  =  1,...,2N)  andDiu=Duk  (k,l  =  l,...,Af). 

In  the  case  of  very  small  excursions  from  the  steady  state  equilibrium  position.  Equations  (6) 
may  be  further  reduced  by  omitting  the  quadratic  nonlinearides  to  give  the  linear  equations 

J  ^ 

MODAL  DECOMPOSITION  OF  LINEAR  EQUATIONS  OF  MOTION 
EIGENVALUES  AND  EIGENVECTORS 

The  natural  frequencies  and  corresponding  mode  shapes  of  the  towed  system  motions 
described  by  the  linear  Equations  (9)  are  determined  by  c^culating  the  2N  eigenvalues  and 
eigenvectors  of  the  matrix  A  whose  elements  are  the  A^^  (ij  =  1,...,2N).  Typically,  some  of  the 
eigenvalues  will  be  real,  and  the  remaining  will  occur  in  complex  conjugate  pairs.  Assiuning  that 
the  eigenvalues  all  have  negative  real  parts  (indicating  that  the  steady  state  equilibrium  is  linearly 
stable),  the  real  eigenvalues  represent  overdamped  modes,  and  the  complex  ones  rejm^nt 
underdUunped  modes.  It  will  be  assumed  in  tiie  balance  of  tius  report  that  the  steady  state  equilibrium 
is  linearly  stable  and  that  A  has  2N  distinct  eigenvalues  and  hence,  2N  linearly  independent 
eigenvectors. 

As  discussed  in  Refnence  15,  the  eigenvector  associated  with  the  eigenvalue  A,  takes  the 

form 


9i 


IWp*")/ 


(10) 


where  is  a  2N  x  1  vector  that  can  be  partitioned  into  two  Nxl  vectors;  as  indicated  in  Equation 
(10),  the  bottom  half  vector  is  simply  the  product  of  the  eigenvalue  and  the  top  half  vector.  The 
eigenvector  is  real-valued  for  real  eigenvalues  and  complex-valued  for  complex  eigenvalues. 
Mtneover,  eigenvectors  that  correspond  to  a  pair  of  complex  conjugate  eigenvalues  are  Aemselves 
conqrlex  conjugates  of  each  other. 


The  top  half  of  the  eigenvector  represents  the  shape  of  the  towed  system.  If  the  eigenvector 
is  teal-vidue^  then  all  points  along  the  cable  move  in-phase  as  the  mode  is  excited.  However,  if 
the  eigenvector  is  complex,  then  there  is  a  phase  shift  between  the  motions  at  different  locations 
along  the  cable.  Although  the  phase  angle  does  vary  along  the  cable,  the  relative  phase  shift  between 
any  two  points  along  the  cable  is  truiintained  as  the  ctdile  rrxrves.  For  example,  it  is  not  uncommon 
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fen*  motion  near  the  system  tow  point  to  be  180  degrees  out  of  phase  with  the  cable  near  the  towed 
body.  This  is  illustrated  in  Figtm  2  that  shows  the  steady-state  configuration  and  five  snap-shots 
of  an  in-plane  mode  of  a  100-foot  cable  towing  a  SOO-pound  spherical  body  at  5  knots.  *^0  dis¬ 
placements  have  been  magnified  for  clarity. 

Finally,  P,  the  "eigencolumn"  or  "modal”  matrix  of  A  is  defined  to  be  the  2^  x  2A^  matrix 
whose  columns  are  fomiffid  from  the  eigenvectors  of  A.  When  the  eigenvalues  are  all  distinct,  the 
modal  matrix  can  be  used  to  diagonalize  A  as  follows: 

r‘AP  =  [A].  (11) 

Here  P  represents  the  nKxlal  matrix  and  [A]  represents  a  diagonal  matrix  that  has  the  eigenvalues 
of  A  on  its  diagonal. 


SOLUTION  OF  LINEAR  EQUATIONS  OF  MOTION 

The  linear  equations  of  motion  (Equation  (9))  can  now  be  solved  by  introducing  a  set  of 
normal  or  modal  coordinates  that  are  related  to  the  state  variables  through  the  modal  matrix  as 
follows: 


y.  =  (i,m  =  l,...,2N)  (12) 

m 

where  the  are  the  elements  of  the  modal  matrix  P.  Substituting  this  relationship  into  the  linear 
equations  with  external  inputs  and  premultiplying  that  equation  by  leads  to  a  set  of  uncoupled 
differential  equations  for  the  modal  coordinates  that  may  be  written  as  follows: 

f,  =  Vr+/r  (r  =  l,...,2N)  (13) 

where 

f.  =  (14) 

*  »,* 

Equation  (13)  represents  a  set  of  uncoupled  complex- valued  first-order  ordinary  differential 
equations  for  the  z,.  As  shown  in  Reference  12,  the  general  solution  to  Equation  (13)  can  be  found 

by  using  the  integrating  factor  the  solutitm  may  then  be  written  in  either  of  the  two  equivalent 

forms 


+  (r  =  l . 2)V) 

(15) 

+  e'^  (' e'*’'f,(x)dx  (r  =  l . 2)V). 

Jo 

(16) 
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Here  the  =  represent  2N  arbitrary  constants  of  integration  and  the 

z,(0)  (r  =  I,..  .,2N)  are  the  values  of  the  z,  at  time  r  =  0.  The  values  of  the  state  variables  can  now 
be  recovered  by  substituting  the  right  sides  of  either  Equation  (IS)  or  (16)  into  Equation  (12)  to 
give 


or 


y.  =  j  e'^y,(x)dx 


yi  = 


t 


(17) 

(18) 


Equations  (17)  or  (18)  give  the  general  solution  to  the  linear  dynamical  equations  of  motion  of  the 
N  degree  of  freedom  finite-segment  model  of  a  towed  system  with  arbitrary  external  excitation. 


ANALYSIS  OF  SYSTEM  MOTIONS  BY  MODE  SUPERPOSITION 

Motions  of  a  towed  system  that  can  be  described  by  the  linear  equations  represented  in 
Equation  (9)  can  be  written  as  linear  combinations  of  the  eigenvectors  of  that  system.  First,  z  (0), 
the  initial  modal  vector  is  determined  from  y  (0)  the  initial  state  vector  as  follows: 

m  =  F'y(0).  (19) 

The  entries  of  the  modal  vector  can  then  be  substituted  into  Equation  (18)  to  determine  the  system’s 
motion  throughout  time. 

Note  that  because  the  motions  represented  by  the  yj  on  the  left  sides  of  Equation  (18)  are 
real-valued,  the  sum  on  the  right  side  must  also  be  real.  As  discussed  in  Reference  12,  this  requires 
that  underd^p^  modes  be  excited  in  complex  conjugate  pairs.  (Corresponding  to  each  such  pair 
of  modes  is  a  pair  of  modal  coordinates  that  are  themselves  complex  conjugates  of  each  other.  So 
then,  if  only  p,  and  Py,  a  single  pair  of  underdamped  (complex  conjugate)  modes  are  excited,  the 
resulting  time-varying  state  vector  can  be  written  as  follows: 

y(r)  =  z.(0pi  +  zy(r)py  =  ZiPj+i^p;  =  ZiPi  +  fei^  (20) 

where  z,  and  Zy  represent  the  two  complex  conjugate  modal  coordinates  associated  with  p,  and  py 
and  the  overbars  represent  complex  conjugates.  Becau%  the  right  side  of  Equation  (20)  is  the  sum 
of  a  complex  conjugate  pair,  the  resulting  state  vector  is  real- valued.  In  general,  therefore,  the  sum 
presented  by  Equation  (18)  will  always  be  real- valued  since  complex  values  always  appear  in 
complex-conjugate  pairs. 
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Note  also  that  for  systems  undergoing  external  excitation,  the  above  observations  indicate  an 
underlying  structure  for  the  B  array  defined  in  Equation  (14).  In  particular,  if  z,  and  Zj  represent  a 

pair  of  complex  conjugate  modal  coordinates,  then  B^  =  {Bj^)  {k  =  to  ensure  proper 

excitation  of  both  modes. 


REDUCED  NONLINEAR  EQUATIONS  IN  MODAL  VARIABLES 


Using  the  techniques  outlined  in  the  previous  section,  the  reduced  nonlinear  Equations  (6) 
can  also  be  transformed  into  modal  variables.  This  is  accoixq>lished  by  substituting  from  Equation 
(12)  into  Equation  (6)  and  premultiplying  by  and  simplifying.  The  result  is  a  set  of  differenaal 
equations  for  the  modal  coordinates  that  may  be  written  as  follows: 


i.  =  U  +  SB*v*+^lC„z^,  +  ^lD*,v*v,+  lE^v*z„  (r  =  l,...,2N)  (21) 

k  Z,T,*  Z  k,l  m,k 


where 


«.  =  (22) 

I 

»  tj 

Equation  (21)  represents  a  set  of  complex-valued  first-order  ordinary  differential  equations 
for  the  z,  that  are  coupled  through  the  quadratic  nonlinearities.  Unlike  the  linear  uncoupled  equations 
represented  by  Equation  (13),  these  equations  are  not  so  easily  solved.  To  better  understand  the 
structure  of  these  equations,  some  comments  on  the  properties  of  the  C,  D,  and  E  arrays  may  be 
helpful  at  this  point.  As  with  the  C  and  D  arrays  of  Equation  (6),  the  C  and  D  arrays  possess  the 
following  symmetry  properties: 

=  {f,r,s  =  1,...,2N)  (r  =  1 . 2N;  k,l  =  1 . Af).  (24) 

Also,  as  a  result  of  the  fact  that  the  modal  variables  associated  with  pairs  of  complex  conjugate 
modes  must  them^lves  be  complex  conjugates  of  each  other,  addition^  symmetries  must  exist  for 
the  B,  (^,  D,  and  E  arrays.  In  particular,  if  and  P;  represent  a  pair  of  complex  conjugate  modes, 
then  the  following  relationships  hold: 

=  (25) 

D^={Piu)  (t./  =  l . M)  (26) 

£a=(£5  (t  =  l . *f).  (27) 

To  fully  understand  the  effects  of  quadratic  nonlinearities  on  the  modal  response  of  underwater 
towed  systems.  Equations  (21)  must  be  solved  as  2N  coupled  equations.  As  widi  the  fully  nonlinear 
Equations  (1)  and  (2),  this  can  be  done  by  numerical  integration;  unfortunately,  anlayticd  solutions 
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to  these  equations  are  generally  not  available.  This  approach  provides  perfectly  acceptable  solu¬ 
tions;  however,  to  gain  significant  insight  into  the  fundamental  dynamical  behavior  of  these  systems, 
many  such  numerical  solutions  must  be  generated  and  studied.  This  trial-and-error  process  is  very 
tedious. 

Actually,  a  close  examination  of  Equations  (21)  for  sample  systems  shows  that  these  equations 
are  not  fully  coupled.  (Consider,  for  example,  the  towed  system  of  Figure  2,  consisting  of  a  sphere 
and  a  100-foot  cable  being  towed  at  5  knots  with  no  side  currents.  The  modes  for  this  system  may 
be  separated  into  two  basic  categories:  1)  modes  that  describe  motions  in  the  vertical  plane  of  the 
mean-ship  motion  (referred  to  here  as  in-plane  modes;  and  2)  modes  that  describe  motions  of  the 
system  out  of  this  plane  (referred  to  here  as  out-of-plane  modes).  Inspection  of  the  numerically 
generated  C  matrix  for  this  case  shows  that  the  coefficients  coupling  some  of  the  in-plane  and 
out-of-plane  nKxles  are  small,  indicating  the  possibility  that  coupling  between  those  modes  can  be 
ignored. 

Unfortunately,  even  taking  advantage  of  these  simplifications  still  leaves  one  with  too  many 
equations  for  analytical  methods.  However,  if  these  equations  are  further  simplified  by  ignoring 
most  or  all  of  the  mode  coupling,  perturbation  (asymptotic)  expansions  may  be  sought  to  provide 
approximate  analytical  solutions.  To  this  end,  the  following  section  describes  how  the  method  of 
multiple  scales  may  be  applied  to  the  equation  for  a  single  underdamped  mode  and  notes  some  of 
the  insights  gained. 


NONLINEAR  PERTURBATION  ANALYSIS 

As  mentioned  above,  the  nonlinear  dynamic  response  of  towed  systems  can  be  modeled  by 
nonlinear  equations  of  the  form  of  Equations  (1)  and  (2).  If  we  neglect  cubic  and  higher  nonlin¬ 
earities,  Equations  ( 1 )  and  (2)  can  be  r^uced  to  the  form  shown  in  Equations  (6),  and  these  equations 
can,  in  turn,  be  transformed  into  a  set  of  equations  for  modal  (or  normal)  variables  as  shown  in 
Equation  (21).  Even  though  these  equations  are  not  fully  coupled,  application  of  analytical  per- 
tui^ation  techniques  to  solve  these  equations  is  impractical.  Hence,  to  apply  procedures  such  as 
the  method  of  multiple  scales  to  find  approximate  solutions  to  the  equations,  furAer  simplifications 
must  be  made.  Obviously,  simplifications  made  without  sound  jus^cation  may  lead  to  erroneous 
conclusions;  consequently,  any  conclusions  drawn  fiom  such  andyses  must  be  verified  by  returning 
to  the  numerical  solution  of  the  fully  nonlinear  Equations  (1)  and  (2). 

As  a  first  step,  it  will  be  assumed  that  all  modes  (except  for  complex-conjugate  pairs)  are 
decoupled,  allowing  a  single  equation  to  be  written  for  each  mode.  The  following  sections  outline 
the  form  of  the  equation  for  an  underdamped  mode  with  tow  point  excitation  and  the  application 
of  the  method  of  multiple  scales  to  develop  an  asymptotic  solution  to  the  equation. 


EQUATION  FOR  AN  UNDERDAMPED  MODE 

Ignoring  the  coupling  between  separate  modes  and  noting  the  symmetries  that  exist  in  the  B , 
C,  15,  and  t  arrays  given  in  equations  (25)  through  (27)  for  modal  variables  corresponding  to 
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complex  conjugate  inodes,  the  equation  for  a  single  underdamped  mode  of  a  towed  system  excited 
by  tow  point  motion  pertuibations  can  be  shown  to  be  of  the  following  form: 

dz  _  A. 

=  8^,  +  <v*+p,5^+Y^,J,+i(t)z,+toz,+Ti,(x)  (28) 

where 

ximrutanensioruil  time)=^\'kr\t  8,=VlXyl  (29) 

a,=((;„y2|Xj  P,=((;„y2ix,i  y=(.C:„+Cj)rAVt=f.C„y\^  m 

where  the  indices  "r  "  and  ”s "  represent  the  complex  conjugate  mode  numbers,  and  overbars  represent 
complex  conjugates.  For  harmonic  tow  point  excitation,  the  forcing  terms  ^^(x).  Vr(x)t  and  Ti,.(T) 
may  be  written  in  the  following  form: 

to  =  to = OD 

to (32) 

where  the  nondimensional  frequency  (o  is  related  to  the  forcing  frequency  q  as  follows:  cd = ^  X,| . 
The  values  of  the  and  the  in  the  above  equations  are  functionally  dependent  on  the  excitation 

amplitude  and  frequency  and  elements  of  the  E  array;  the  are  functionally  dependent  on  the 
excitation  amplitude  and  frequency  and  elements  of  the  B  and  D  arrays.  (Consequently,  the 
the  and  the  i$„(r  =  1, 2)  are  coefficients  of  terms  linear  in  the  external  input  pertuibations,  and 

the  are  coefficients  of  terms  that  are  quadratic  in  the  external  input  perturbations. 

MULTIPLE  SCALES  ANALYSIS 

To  determine  an  an  approximate  solutitm  to  Equation  (28),  the  method  of  multiple  scales‘s 
is  employed.  To  this  end,  a  solution  is  sought  for  sn^  but  finite  amplitudes  of  the  form 

Zr  =  ^rl  +e’^rJ  (33) 

with  the  time  scales  Tj=efx(i  =:0,1,2).  To  maintain  the  terms  of  ti,.(x)  that  are  linear  in  the  external 
input  perturbations  as  first-order  effects  in  this  expansion  uid  all  other  terms  as  second-order  effects, 
it  will  be  assumed  that  the  forcing  terms  $r(x)«  Vr(x),  and  tj^x)  may  be  written  in  the  following  form; 

to=e(Gw«‘**+Qr2«‘'*‘)=o>,(x)  to=e(ew«‘*'+G.2«-"*)=eVr('')  (34) 
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Substituting  from  Equations  (33)  through  (35)  inro  the  differential  Equation  (28)  and  separating 
equations  based  on  like  powers  of  £  gives  the  following  three  partial  differential  equations  for  the 
z^(i*  1,2,3): 

(36) 


+P^;  +*rCr.)Zrf+ +’wcr(i)  (?■) 

dz,j  dz,2  3z,;  _  _  _ 

^-8^,2  +  +  Y,(z,/2,2  +  ZrfZrz)  + 


<|),(ro)2,2  +  V,(7’o)^,2.  (38) 


Using  Equation  (16),  the  solution  to  Equation  (36)  can  be  written  as  follows: 

where 


A=A{T^,T^ 


B  = 


-Rrl 

8r-it0 


C  = 


-Rr2 

6^+ico' 


(39) 

(40) 


As  usual  with  the  method  of  multiple  scales,  the  coefficient  A  of  the  homogeneous  solution  is 
assumed  to  be  a  function  of  the  slower  time  scales.  The  dependence  of  A  on  these  scales  will  be 
determined  at  subsequent  levels  of  the  perturbation  expansion  through  the  elimination  of  secular 
terms. 

Substituting  now  from  Equation  (39)  into  the  right  side  of  Equation  (37)  and  using  the  second 
term  m  the  right  side  of  Equation  (16)  yields  the  following  equation  for  the  particular  solution  to 
Equation  (37): 
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where  the  coefficients  Dy  through  Dyy  may  be  written  as  follows: 


D,=<v4*  D^=yM  (42) 

D^=2a^B+y^C+AQ,,  Di=2xi^C +y^B+AQ^2  (43) 

Z)«=2P;4C+y;4B+JQ,7  Z),=2p;4B+Y;4C+Je^  (44) 


D,=ovB*+ Y^C  +BQ,,  +  CQ„  +R,3  (45) 

D,  *  a,C*+  p^+ Y^C + Cdi  +BQa  +Rr3  (46) 

D,o = 2a^C + 2pJC+y,(BB  +  CC) + CQ,/  +5^2  +^2./ + (47) 


-Dy  -D2 


— D3  — D4  D5  -:Dj  — D7 

2S,-8,  B,  ito  ica  ^-5,+i© 


^9  ^10 

6,-i2m  8,.+i2a)  5,  ’ 


(48) 


Eliminating  the  secular  term  firom  this  equation  requires  that  (dAf^Ty) = 0  and  consequently  that 
A  =ACrd.  It  should  also  be  noted  here,  however,  that  this  solution  is  singular  whenever  any  of  the 
denominators  in  Equation  (41)  become  zero;  this  occurs  when  ©=23(5,.)  which  makes 

-  5,. + 1© = 0.  This  condition  corresponds  to  the  condition  that  © = 23(X,.),  that  is  that  the  forcing 
ftequen^  is  twice  the  ficequency  of  finee-modal  response.  An  expansion  valid  in  the  neighborhood 
of  this  singularity  is  derived  in  the  following  section. 
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To  determine  how  A  depends  on  Tj  and  to  ensure  an  asymptotic  expansion  valid  through 

0(e^),  one  needs  only  eliminate  the  secular  terms  from  the  right  side  of  Equation  (38).  Substituting 
fnm  Equation  (39)  and  Equation  (41)  into  Equation  (38)  and  setting  the  sum  of  the  secular  terms 
to  zero  gives  the  following  equation  for  the  coefficient  A : 

(49) 


where 


^  f-D^o  BE,  C£4)  .of  BEj  CE^  BE. 

(  ^£7  ^  BE4  CEj  ^  -^rlBs  ^  Qr2^4  ^  QtJ^b  ^ 

\  5, -5, +1(0  i(o  1(0  ^  ICO  ICO  8^ -^-1(0  ^5,-^+ 


Ir 


£4s2oc^+Y,C+ai  £5=2(X,C+%£+G,2 

£j*2p,C  +7^  +Q,/  £,=2p^  +Y,C  +  Q^. 

Therefore,  an  asymptotic  expansion  for  z,  can  be  written  as  follows: 


2br-br  8^  ICO  1(0  ^ 

tT  ^6  (®r'*’**)ro  ^7  pr~**y*  iJaT,  ~D, 

^8t.-8,+ioo^  ^8,-8i.-i(o*  ^8,.-i2co^  ^8,+i: 


a 


.28/, .  Di  a/, .  i 
^  +■"5 — +■ 

28,-8^ 


_  •2«Dr,  ~D9  - 

i2co^  ^8,  +  i2(0* 


As  with  the  particular  solution  shown  in  Equation  (41),  this  solution  is  singular  when  co  =  23(8,.). 


For  values  of  co  not  in  the  neighlxnliood  of  23(8,.),  Equation  (SO)  can  be  used  to  determine 
numerical  vali»s  for  a,  and  hence,  whether  the  perturbation  expansion  is  bounded  or  unbounded. 
For  example,  should  the  real  part  of  a  be  positive  for  some  range  of  excitation  amplitudes  and 
frequencies,  it  indicates  that  the  free-oscillation  terms  are  large.  There  are  two  possible  explanations 
for  this  occurrance.  The  first  is  that  the  expansion  has  simply  not  been  carried  out  to  a  high  enough 
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order  to  provide  a  uniform  expansion.  The  second  is  that  the  solution  for  z,i  in  Equation  (39)  is 
fundamentally  different  in  character  from  the  exact  solution  to  Equation  (28),  indicating  that  the 
solution  itself  may  exhibit  relatively  strong  nonlinear  characteristics. 

Finally,  note  that  Equation  (S3)  contains  the  constant  term  -D,oE^/5,  that  is  indicative  of  a 
nonsymmetric  cable  response  about  its  steady-state  configuration.  This  behavior  was  noted  in 
References  1  and  S;  it  is  a  direct  result  of  the  quadratic  nature  of  the  fluid  drag  forces  on  the  system. 
The  fluid  forces  on  the  system  increase  faster  as  the  system  moves  in  the  diction  of  the  tow  than 
they  decrease  as  it  moves  opposite  to  the  direction  of  the  tow.  D,o  has  been  found  to  be  nonzero 
for  in-plane  modes  and  zero  for  out-of-plane  modes  as  expected. 


MULTIPLE  SCALES  ANALYSIS  -  EXPANSION  NEAR  SINGULARITY 

As  mentioned  above,  the  expansion  shown  in  Equation  (S3)  is  not  valid  in  the  neighbortiood 
of  (0  =  23(5,).  To  develop  an  expansion  that  is  valid,  introduce  a  so-called  detuning  parameter  a 

defined  such  that  (a=i(B,-S,)-i-eo=23(S,)-l-e(T.  Using  this  definition  converts  the  term  on  the 
right  side  of  Equation  (37)  responsible  for  the  singularity  into  a  secular  term.  In  particular, 

(54) 


Then,  elimination  of  secular  terms  leads  to  the  equation: 

=  0.  (55) 

Assuming  that  the  secular  terms  have  been  omitted,  an  asympmtic  expansion  for  z,  through  0(e^) 
becomes 


z,  =  (£4  (7,,  Tj) + e*Z)„)e*'^*  +  eBe'"*’  -h  eCe 


•«r. 


i 


“*■*0  .  Ji  ~^»0 

'  HT, 


6,  25,-5^  br  <©  to 


D^  -D,  iMT,  -Df 


(56) 


The  dq)endence  of  A  on  the  time  scales  Tj  and  Tj  is  determined  by  the  elimination  of  all  secular 
terms.  As  before,  the  elimination  of  secular  terms  from  the  right  si^  of  Equation  (38)  leads  to  the 
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equation 


dA 

dTi 


(57) 


where  a  is  now  defined  to  be 


BEs  CE^ 

]+  +% 

-tDio 

CEn 

BE^ 

CEs" 

1  8^ 

100  1(D 

> 

-5, +100 

^  10) 

'lOO^ 

f  BEn  ' 


Qrl^S  .  .  Qil^l 

•¥  :  -T* 


1(0  1(0  6i.-^+i(o' 

Together  Equations  (55)  and  (57)  can  be  shown  to  be  equivalent  to  the  differential  equation 

~(eE^‘^)A-(ea)A=Q. 


(58) 


(59) 


This  equation  is  a  first-order  equation  with  tune-varying  coefficients;  to  transform  it  to  an  equation 
with  constant  coefficients,  let  A  Making  this  substitution,  separating  into  real  and 

imaginary  parts,  and  solving  the  coupled  first-order  equations,  give  the  following  result: 


91(G) 


3(G)  = 


'M2-E*9l(fl)-E9l(£:«)' 

f-e«!(a)+c3(£a  ^ 

gi*  + 

f-e?3(a)+e3(£a  ^ 

IV 


(60) 

(61) 


where 


H,2=e*9l(a)±eVp  if  p=|£y*-o*/4+(ea+^)3(a)^0 
m2*e*9l(a)±eiV^  if  p<0. 

The  first  free  oscillation  term  in  this  expansion  grows  without  bound  as  x  ->  <»  whenever  the 
amplitude  A  is  unbounded.  This  (xxnirs  when  p>0  and  ^(a)±Vp>0  and  when 
p<0  and  Sl(a)>0.  As  before,  this  may  indicate  either  that  the  perturbation  solution  has  not 
been  carried  out  to  a  high  enough  order  to  provide  a  unifcnm  expansion  or  that  the  solution  for  z,i 
in  Equation  (39)  is  fuiulainentally  different  in  character  fiom  the  exact  solution  to  Equatitm  (28). 
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NUMERICAL  EXAMPLES 

As  stated  in  the  previous  secdtm,  one  possible  indication  of  relatively  strong  nonlinear 
response  of  individual  modes  of  the  towed  system  is  the  presence  of  unbounded  free  oscillation 
terms  in  tte  poturbadon  expansions  derived  therin.  In  this  secdon,  numerical  results  are  presented 
for  the  individual  response  of  two  noodes  of  a  towed  sphere  subject  to  tow  point  excitadons.  Results 
are  presented  first  fcv  an  in-plane  mode  of  the  system  and  Aen  for  an  out-of-plane  mode  of  the 
system.  In  pardcular,  values  of  the  arguments  of  the  exponential  funcdon  that  appear  in  the  derinidon 
of  the  coefficient  A  ofthefree-responsetermarepresentedoverthefrequencyrangefromO.Ol  ->  S.O 
radians/secmuL  This  range  is  within  the  respmise  range  of  many  surface  ships  and  includes  the  free 
refuse  frequency  of  the  selected  modes. 

The  to\^  system  consists  of  a  IS-inch  diameter  sphere  towed  at  IS  knots  at  the  end  of  ISO 
feet  of  a  O.S4-inch  diameter  smooth  jacketed  cable.  The  sphere  weighs  SOO  pounds  in  air;  thecable 
weighs  0.39  pounds/foot  in  air  and  has  a  buoyancy  of 0.089  pounds/foot  in  sea  water.  For  modeling 
purposes,  the  cable  was  divided  into  10  linlu  (finite-segments)  with  the  longest  link  at  the  towed 
end  of  the  cable  and  the  shortest  link  at  the  sphere;  this  is  consistent  with  the  mode  shapes  of  the 
cable. 


The  lowest-fiequency  underdamped  mode  of  this  system  is  an  out-of-plane  mode  and  the 
next  lowest-fiequency  underdamped  mode  is  an  in-plane  mode.  The  eigenvalues  Xy  associated  with 
these  modes  are 


Out-of-Plane  Mode;  \  s  (-0.69357;  1.2480) 

In-Plane  Mode:  Xy « (-1.0540;  1.4405). 

The  corresponding  mode  shapes  are  animated  in  Figure  3  for  one  complete  cycle  of  free  response 
assuming  no  ampUtude  decay.  For  the  out-of-plane  mode,  the  motion  is  distributed  over  die  entire 
lengdi  of  cable;  larger  amplitudes  are  experienced  at  the  free  end.  In  the  case  of  the  in-plane  mode, 
most  of  the  motion  occurs  at  the  free  end  of  die  cable;  little  motion  occurs  near  the  tow  point 

For  this  system,  the  reduced  nonlinear  ^nations  (6)  were  generated,  transformed  into  the 
fonn  of  Equation  (21),  and  then  reduced  to  a  single  equation  of  the  form  of  ^nation  (28)  for  each 
mode  and  excitation  type.  The  external  excitation  was  applied  both  in  and  normal  to  the  towing 
plane,  the  former  being  referred  to  as  in-pUme  excitation  ^  the  lattn^  out-of-plane  excitation.  The 
in-plane  excitatiem  was  further  subdivided  into  excitations  along  and  noin^  to  the  cable.  The 
non-zero  coefficients  in  Equation  (28)  for  each  mode  and  excitation  are  summarued  in  the  following 
table. 
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NON-ZERO  O 

OEFFICIENTS 

MODE 

In-Plane 

Excitation 

Out-of-Plane 

Excitation 

In-Plane 

a..P,.Yr.«r 

ar.pMYr.'Hrr 

Out-of-Plane 

A  A 

A 

T)rJ 

Note  that  in  the  case  of  out-of-plane  excitation  of  the  out-of-plane  mode,  no  quadratic  nonlineaiities 
are  present;  hence,  the  resulting  equation  is  linear.  Ftx-  in-plane  excitation  of  this  mode,  quadratic 
nonUnearities  arise  only  throu^  the  forcing  terms.  The  equations  for  the  in-plane  mode  are  more 
complex.  For  in-plane  excitation  of  the  in-plane  noode,  all  pmsible  quadratic  nonlinearides  are 
present  For  out-of-plane  excitation  of  this  mcxk,  the  linear  forcing  terms  and  some  of  the  nonlinear 
rarcing  terms  are  absent 

Consider  first  results  for  excitation  of  the  in-plane  mode.  Figure  4  shows  the  real  part  and 
Figure  5  shows  the  imaginary  part  of  the  exptMient  a  the  expansion  in  Equation  (S3)  for  excitatitm 
pandlel  to  the  towed  end  of  tte  cable;  the  excitation  amplitude  is  0.5  feet  (1  foot  peak  to  trough) 
and  is  applied  over  the  frequency  range  0.01  ^  m  ^  5.0  radians/second.  Note  that  results  for  the 
range  2.8  <!  m  ^  3.0radians/second  have  been  mnitred  to  avoid  the  singularity  at  m = 2S(X,.)  =  2.881 
radians/second.  Even  though  the  real  part  of  a  appears  unaffected  by  the  presence  of  the  singularity, 
it  is  clear  from  Figure  5  that  the  imaginary  part  of  a  is  affected  by  the  presence  of  the  singularity. 
In  fact,  if  one  plots  the  teal  and  imaginary  parts  of  a  over  the  range  2.8  ^  ^  3.0  radians/second, 
it  is  clear  that  both  values  are  affect^ 

For  an  expansion  in  the  neighborhood  of  this  frequency,  an  expansion  of  the  form  of  Equation 
(56)  is  sought;  the  coefficient  >4  must  satisfy  Equation  (59).  Unforumately,  as  shown  in  Figure  6, 
the  teal  part  of  is  positive  (albeit  small)  over  this  range.  However,  as  fte  excitation  amplitude 
is  lower^  the  frequency  range  over  which  the  exponents  have  positive  real  parts  becomes  smaller 
and  strudler  around  tiie  singitiarity.  In  conclusitm,  the  perturbatitm  expansions  devel(^)ed  above 
m  branded  throughout  the  frequency  range  shown,  except  in  the  neighbc^ood  of 
ms23(X^)- 2.881  radians/second.  BouiKled  solutions  can  be  found  closer  and  closer  to  this 
frequency  by  lowering  the  excitation  amplitude. 
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Next,  consider  excitation  of  the  in-plane  nsode  perpendicular  to  the  towed  end  of  the  cable. 
Figure  7  shows  the  real  part,  and  Figure  8  shows  the  imaginary  part  of  the  exponent  a  of  the  expansion 
in  Equation  (S3)  for  an  excitation  amplitude  of  0.5  feet  (1  foot  peak  to  trough)  over  the  fr^uency 
range  0.01  m  S.O  radians/secons.  Note  that,  as  alx)ve,  results  for  the  range  2.8  ^  oo  ^  3.0 
radians/secons  have  been  omitted  to  avoid  the  singularity  when  m  =  22(X)  =  2-88 1  radians/second. 

Unlike  the  above  case,  the  real  part  of  a  is  positive  over  the  entire  range  of  Oequencies  shown. 
Moreover,  lowering  the  excitation  amplitude  to  0.001  feet  simply  lowers  the  values  of  the  exponents, 
but  they  remain  p^tive.  Hence,  even  for  extremely  small  excitation,  the  perturbation  expansion 
is  unbonded,  llus  is  true  in  the  neighborhood  of  the  singularity  as  well. 

The  results  for  out-of-plane  excitation  of  this  mode  are  shown  in  Figure  9.  The  real  part  of 
a  is  negative  over  the  entire  frequency  range,  indicating  a  bounded  expansion.  Note  that,  unlike 
the  case  of  in-plane  excitation,  no  singularity  exists  in  this  case  because  the  coefficient^^  in  Equation 
(41)  is  zero,  hence  eliminating  the  singular  term  in  the  expansion. 

Consider  finally,  results  for  excitation  of  the  out-of-plane  mode.  Since  Equation  (28)  reduces 
to  a  linear  equation  for  out-of-plane  excitation  of  the  out-of-plane  nxxle,  pntivbation  methods  are 
not  required  to  study  its  response  when  cubic  and  higho’  order  nonlinearides  are  ignored.  For 
in-plaire  excitation  of  this  mode,  the  steady-state  solution  is  zero  and  nonzero  response  exists  only 
for  iKMizero  initial  conditions;  that  is,  the  noodal  coordinate  simply  relaxes  from  its  perturbed  initi^ 
state  back  to  zero.  For  in-plane  excitation  along  or  normal  to  the  cable,  the  real  part  of  a  is  zero 
for  all  frequencies,  and  the  imaginary  part  of  a  is  near  zero  except  very  near  the  singularity  at 
Q  -  23(Xy).  Hence,  the  perturbation  expansion  is  bounded  at  all  frequencies,  except  at  fiequencies 

very  near  the  singularity.  A  plot  of  the  imaginary  part  of  a  for  in-plane  excitation  normal  to  the 
cable  is  shown  in  Figure  10.  As  in  the  cases  atove,  using  detuning  around  the  singularity  still 
resulted  in  an  unbounded  expansion.  Also,  as  before,  as  the  aixq>litude  of  the  external  excitation  is 
lowered,  bounded  expansions  can  be  found  closer  to  the  the  singularity  point 

As  a  cross-check  of  the  above  results,  the  differential  equations  governing  the  noodal  variables 
were  numerically  integrated  at  selected  excitation  frequencies.  The  results  of  the  integration  of  the 
equaticms  (of  the  form  of  Equation  (28))  ctmtaining  the  quadratic  nonlinearities  were  compared 
with  those  of  the  integration  of  the  linetuized  forms  of  these  equations.  These  comparisons  led  to 
the  following  conclusions.  First,  in  cases  where  9t(a)  is  negative,  the  differences  between  the  linear 
and  nonlinear  solutions  were  relatively  small  in  comparison  to  the  cases  when  9t(a)  is  positive. 
Also,  as  9t(a)  gets  larger  (more  positive),  the  nonlinear  behavior  becomes  stronger.  In  gen^,  the 
results  suggest  that  SH(a)  can  be  used  to  indicate  how  nonlinear  the  response  of  individual  nxxles 
are  to  harironic  excitation.  This  is  consistent  with  the  notion  that  perWbation  methods  may  be 
applied  most  successfully  to  weakly  nmdinear  systems. 

It  should  be  no^  here  that  numerical  integration  uncovered  no  particularly  unusual  behavior 
at  the  singular  point  m  =  23(Xy).  As  a  result,  this  singularity  should  not  be  considered  an  indication 
of  highly  nonlinear  behavior  but  rather  a  product  of  the  perturbation  expansimi  itself. 
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SUMMARY 

This  paper  presents  the  methodology  for  developing  perturbation  expansions  to  describe 
the  response  of  individual  nK)des  of  a  towed  system  to  tow  point  excitation  using  the  method  of 
muldpte  scales.  Hie  equations  of  perturbed  motion  (up  through  quadratic  nonlinearities)  are 
develc^ied  about  a  general  steady-state  equilibrium  configuration  using  a  numerical,  finite- 
segnoent  computer  program  used  at  the  Qiastal  Systems  Station  to  model  die  nonlinear  dynamics 
of  towed  cable  systems. 

Results  are  shown  for  two  modes  of  a  sphere  towed  at  IS  knots  from  a  ISO-foot  cable;  one 
of  the  modes  is  characterized  by  motions  in  die  towing  plane  (in-plane)  and  the  other  is  charac¬ 
terized  by  motions  normal  to  the  towing  plane  (out-of-plane).  Each  m^  was  subjected  to  in¬ 
plane  and  out-of-plane  harmonic  excitation  at  the  system’s  tow  point  The  results  indicate  that 
the  out-of-plane  nKxle  responds  Imearly  to  out-of-plane  excitation,  but  that  the  in-plane  mode  has 
non^ear  response  to  in-plane  excitation.  For  the  towed  sphere  system  presented,  the  strength  of 
the  nonlinear  response  depends  on  the  direction  of  die  in-plane  excitation;  excitation  normal  to 
the  towed  end  of  the  cable  is  shown  to  produce  a  stronger  nonlinear  response  than  excitation  par¬ 
allel  to  the  cable.  Weak  nonlinear  lespmise  is  charactnized  by  a  bounded  perturbation  expan- 
si(Mi  and  strong  nonlinear  response  is  characterized  by  an  unbounded  expansion. 

Perturbation  expansions  for  in-plane  excitation  of  in-plane  modes  are  shown  to  have  a  sin¬ 
gularity  when  the  excitation  fiequency  is  twice  die  fiee-res^nse  fiequency  of  that  mode.  The 
fiec^uency  range  over  which  the  expansion  breaks  down  can  be  narrowed  by  lowering  the 
exatation  amplitude. 

Numerical  integration  of  the  differential  equations  governing  the  modal  variables  (associ¬ 
ated  with  die  above  modes)  supports  the  ctmclusions  drawn  fiom  Ae  penurbation  analysis. 
However,  no  unusual  behavior  was  found  at  die  frequency  associated  with  the  singularity  in  the 
perturbation  expansion,  indicating  that  the  singularity  was  a  product  of  the  expansion  alone  and 
not  an  indicatin’  of  highly  nonlinw  behavior. 

The  present  analysis  could  be  easily  extended  to  include  systems  with  extonal  excitations 
at  locations  other  than  the  system  tow  point,  such  as,  say,  the  control  surfaces  of  a  towed  vehicle. 
It  can  also  be  extended  to  study  the  coupled  response  of  pairs  of  modes  of  towed  systems  and  the 
effects  of  cubic  nonlinearities.  The  procedure  can  be  us^  to  investigate  the  fundamental 
dynamic  properties  of  towed  systems  and  could  possibly  be  used  to  generate  low-order  dynamic 
inodels  tl^  inclu^  a  minimal  number  of  modes  to  des^be  system  d3mamics.  Numoic^  inte¬ 
gration  of  tile  equations  of  motion  such  models  can  be  compared  with  numerical  integration  of 
the  fully  nortiinear  Equations  (1)  and  (2)  to  determine  tiie  accuracy  of  these  models. 
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FIGURE  1.  TOWED  SYSTEM  WITH  A  SINGLE  VEHICLE  AND  TOW  CABLE 
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FIGURE  2.  MODE  SHAPE  OF  100 -FOOT  CABLE  WITH  SPHERE 

AT  5  KNOTS 
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a.  IN-PLANE  MODE 
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b.  OLT-OF-PLANE  MODE 

nCURES.  MODE  SHAPE  OF  ISO -FOOT  CABLE  WITH  SPHERE 

AT  15  KNOTS 
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RGURE  4.  REAL  PART  OF  EXPONENT  a  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  ALONG  CABLE 
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FIGURE  5.  IMAGINARY  PART  OF  EXPONENT  a  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  ALONG  CABLE 
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FIGURE  6.  REAL  PARTS  OF  EXPONENTS  AND  M2  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  ALONG  CABLE 
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FIGURE  7.  REAL  PART  OF  EXPONENT  a  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  NORMAL  TO  CABLE 
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FIGURE  8.  IMAGINARY  PART  OF  EXPONENT  a  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  NORMAL  TO  CABLE 
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FIGURE  9.  EXPONENT  a  FROM  PERTURBATION  ANALYSIS: 
OUT-OF-PLANE  EXCITATION 
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FIGURE  10.  IMAGINARY  PART  OF  EXPONENT  a  FROM  PERTURBATION 
ANALYSIS:  IN-PLANE  EXCITATION,  NORMAL  TO  CABLE 
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